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In this work, we explore the cosmological consequences of the recently released Constitution 
sample of 397 Type la supernovae (SNIa). By revisiting the Chevallier-Polarski-Linder (CPL) pa- 
rameterization, we find that, for fitting the Constitution set alone, the behavior of dark energy (DE) 
significantly deviate from the cosmological constant A, where the equation of state (EOS) w and 
the energy density pA of DE will rapidly decrease along with the increase of redshift z- Inspired by 
this clue, we separate the redshifts into different bins, and discuss the models of a constant w or a 
constant pa in each bin, respectively. It is found that for fitting the Constitution set alone, w and p^ 
will also rapidly decrease along with the increase of z, which is consistent with the result of CPL 
model. Moreover, a step function model in which pA rapidly decreases at redshift z ~ 0.331 presents 
a significant improvement {Ax^ - -4.361) over the CPL parameterization, and performs better than 
other DE models. We also plot the error bars of DE density of this model, and find that this model 
deviates from the cosmological constant A at 68.3% confidence level (CL); this may arise from some 
biasing systematic errors in the handling of SNIa data, or more interestingly from the nature of DE 
itself. In addition, for models with same number of redshift bins, a piecewise constant pa model 
always performs better than a piecewise constant w model; this shows the advantage of using pa, 
instead of w, to probe the variation of DE. 
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I. INTRODUCTION 



Although it has been a decade since the discovery of the cosmic acceleration |lli|2|], the nature of dark 
energy (DE) still remains a mystery. The most obvious theoretical candidate of DE is the cosmological 
constant A, which can fit observations well, but is plagued with the fine-tuning problem and the coincidence 
problem yl]. Numerous other dynamical DE models have also been proposed in the literature, such as 

1 ^, ta chyon |7], holographic |8], agegraphic |9], hessence I lol . 



quintessence 14)], phantom [5], ^-essence 



Chaplygin gas |11], Yang-Mills condensate 11211 . ect. 

A most powerful probe of DE is Type la supernovae (SNIa), which can be used as cosmological standard 
candles to measure directly the expansion history of the universe. A large sample of nearby Snia with 
z < 0.08 has recently been published il3n . Adding this to the Union sample 11411 leads to the so-called 



Constitution set jlSll which is currently the largest SnIa sample to date. By analyzing this Constitution set 



with Chevallier-Polarski-Linder (CPL) ansatz il6ll . Shafieloo et al. [17] argue that cosmic acceleration may 



have already peaked and that we ai^e currently witnessing its slowing down. By separating the redshifts 
into different bins and assuming a constant equation of state (EOS) w in each bin, Qi et al. pist] also point 
out that this Constitution set shows a deviation at 68.3% confidence level (CL) from the ACDM model. 



Besides, by studying the constraints of DE on the w -w' plane, Chen et al. 11911 claim that ACDM model is 
disfavored by the current observational data at 68.3% confidence. 

Although constraining EOS w of DE is a popular and widely-used method to investigate dark energy, 
Wang and Freese I20I1 pointed out that DE density pa can be constrained more tightly than EOS w given 
the same observational data. Therefore, it would also be very interesting to consider the cases of piecewise 
constant p^. In this work we separate the redshifts into different bins, and discuss the models of a constant 



w or a constant pA in each bin, respectively. Binned fits of EOS have been applied before in 



I 23h . and the Union paper itself [ 14]. Similar analyses have been performed for density binning 12011 . Our 
paper differs from previous ones in two aspects: First, in previous papers the redshift bins are determined by 
setting the discontinuity points of redshift by hand, while in our paper we treat the discontinuity points of 
redshift as free parameters, and find that much smaller ^^^^.^ can be obtained. Second, not only the piecewise 
constant w models but also the piecewise constant pA models are considered in this work. Moreover, by 
comparing these two classes of model, it is found that for same number of redshift bins, a piecewise constant 
Pa model always performs better than a piecewise constant w model. 

This work is organized as follows. In Section 2, we briefly review the models considered here and 
the method of data analysis. In Section 3, we introduce the observational data and describe how they are 
included in our analysis. Section 4 is divided into three parts: First of all, we revisit the CPL model and 
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find that for fitting the Constitution set alone the behavior of DE significantly deviate from the cosmological 
constant A, where w and pA will rapidly decrease along with the increase of redshift z- This result implies 
that the behavior of dark energy might be very different in different slices of redshifts, and inspires us to 
separate redshifts into several bins and to investigate the piecewise constant w model and the piecewise 
constant pA model. Then, we discuss the models with two bins. It is found that for fitting the Constitution 
set alone, pA will also rapidly decrease along with the increase of z, which is consistent with the result of 
CPL model. Moreover, a step function model in which pA rapidly decreases at redshift z ~ 0.331 presents 
a significant improvement (A^^ = -4.361) over the CPL parameterization, and performs better than other 
DE models. We also plot the error bars of pA of this model, and find that this model deviates from the 
cosmological constant A at 68.3% confidence level. Next, we briefly discuss the models with three bins, 
and find that the cases of three bins are very similar with that of two bins. At last, we give a short summary 
in Section 5. In this work, we assume today's scale factor ao = 1, so the redshift z satisfies z = - 1; the 
subscript "0" always indicates the present value of the corresponding quantity, and the unit with c = h = I 
is used. 



II. MODEL AND METHODOLOGY 



Standard candles impose constraints on cosmological parameters essentially through a comparison of 
the luminosity distance from observation with that from theoretical models. In a spatially flat Friedmann- 
Robertson- Walker (FRW) universe (the assumption of flatness is motivated by the inflation scenario), the 
luminosity distance di is given by 

^0 Jo E{z') 

with 

E{z) = H{z)/Ho - [n^oil + z? + (1 - n^o)mf'^ , (2) 

where H{z) is the Hubble parameter, Hq is the Hubble constant, Q^o is the present fractional matter density, 
and f{z) = Pa(z)/pao is a key function, because DE parameterization schemes enter through f{z). 
We will consider the familiar CPL ansatz 1 16], in which the EoS of DE is parameterized as 

z 

w = wo + wi- , (3) 

1 +z 

where wq and wi are constants. As is well known, the corresponding f{z) is given by 

x3(l+W0+Wl) pvr, / ^^^^ 



f(z) = (l+ z)^^'+»'o+-i^ exp ( . (4) 
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For the case where EOS w is piecewise constant in redshift, f(z) can be written as i22l 12311 

n-l 



fiZn-l < Z < Zn) = (1 + {{{l + zO'^'''""''"^'^ (5) 

1=0 

where w, is the EOS parameter in the redshift bin defined by an upper boundary at Zi- This class of 
models has been extensively studied in the literature. As mentioned above, it is interesting to consider the 
piecewise constant pA models. For this case, /(z) can be written as 

f 1 < z < zi 

fiz) - . (6) 

[ e„ Zn-\ <z<z„{n> 1) 

Here e„ is a a piecewise constant, and from the relation ^(O) = 1 one can easily obtain ei = 1. So for 
same number of redshift bins, the number of free parameters of piecewise constant p\ model is one fewer 
than that of piecewise constant w model. It should be mentioned that there are different opinions in the 



2311 . the authors directly set 



literature about the optimal choice of redshift bins in constraining DE. In 11221 
the discontinuity points of redshift as zi = 0.2, Z2 - 0.5, and zs = 1.8. In 12411 . Wang argue that one should 
choose a constant Az for redshift slices. In this work, we just treat the discontinuity points of redshift as 
model parameters in performing the best-fit analysis. As seen below, this leads to a much smaller 

In this work we adopt statistic to estimate parameters. For a physical quantity ^ with experimentally 
measured value ^o, standard deviation cr^, and theoretically predicted value the;^'^ value is given by 

,2 ; 



The total ;t' is the sum of all;^fp, i.e. 



= Yj^I (8) 



Assuming the measurement errors be Gaussian, the likelihood function is given by 

£ oc (9) 



For comparing different models, a statistical variable must be chosen. The ^^.^ is the simplest one, but it 
has difficulty to compare different models with different number of parameters. In this work, we will use 
Xmin^^^f ^ niodel selection criterion, where dof is the degree of freedom defined as 

dof = N- k, (10) 

here N is the number of data, and k is the number of free parameters. Besides, to compare different models 



with different number of parameters, people often use the Bayesian information criterion [|25[1 given by 12611 



BlC = -l\nZ^ax + k\nN, (11) 
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where JLmax is the maxunum likelihood. It is clear that a model favored by the observations should give 
smaller xl,J do f and BIC. 

III. OBSERVATIONAL DATA 
A. Type la supernovae 

For SNja data, we use the latest 397 Constitution sample, the distance modulus fXobsiZi), compiled in 
Table 1 of [ 151 . The theoretical distance modulus is defined as 

l^thizi) = 5 log 10 DtiZi) + Ho, (12) 

where yUo = 42.38 - 5 logjo h with h the Hubble constant Hq in units of 100 km/s/Mpc, and 

is the Hubble-free luminosity distance HqcIl in a spatially flat FRW universe, and here 9 denotes the model 
parameters. The x'^ for the SNIa data is 

397 2 

XsNy^> - 2_t 'Z2. ' ^^^> 



where HobsiZi) and cr, are the observed value and the corresponding Icr error of distance modulus for each 
supernova, respectively. The parameter //q is a nuisance parameter but it is independent of the data and the 



dataset. Following 12711 . the minimization with respect to /iq can be made trivial by expanding the;if^ of 



Eq.(fT4l) with respect to //q as 



xlf,(e) = A(9) - liioBiQ) + i4c, (15) 



where 

iMohsizi) -Hthizn^Q = 0,9)Y 



it 

397 

g^g^ ^ y ^ipbsizd-fithizhfiQ ^0,6/) ^^^^ 

;=1 1 

397 

c 



Evidently, Eq. ([T4l ) has a minimum for - B/C at 



i=i '^i 



xj^(9)=A{9)-^. (19) 

Since = ^5^„„„, instead minimizing one can minimize which is independent of the nui- 

sance parameter fig. 
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B. Baryon Acoustic Oscillations 



For BAO data, we consider the parameter A from the measurement of the BAO peak in the distribution 



of SDSS luminous red galaxies, which is defined as i28h 



where zt = 0.35. The SDSS BAO measurement 



if 

Zh Jo 



dz' 



E{z') 



2/3 



m gives Aohs = 0.469 (W0.98) 



-0.35 



(20) 



0.017, where the 



2911 . It is widely believed that A 



scalar spectral index is taken to be = 0.960 as measured by WMAP5 
is nearly model-independent and can provide robust constraint as complement to SNIa data. The;^;^^ for the 
BAO data is 

{A-Aobsf 



2 

Xbao 



where the corresponding Icr errors is cta = 0.017. 



(21) 



C. cosmic microwave background 



For CMB data, we use the CMB shift parameter R, which is given by 



inO 



dz' 

Eiz'Y 



(22) 



where the redshift of recombination Zrec - 1090 [29]. The shift parameter R relates the angular diameter 
distance to the last scattering surface, the comoving size of the sound horizon at z^. and the angular scale 
of the first acoustic peak in CMB power spectrum of temperature fluctuations |3^ 3]J]. The measured value 



of R has been updated to be Robs ^ 1-710 ± 0.019 from WMAP5 [29]. It should be noted that, different 
from the SNla and BAO data, the R parameter can provide the information about the universe at very high 
redshift. The for the CMB data is 

(R-Robs)^ 



2 

XcMB 



where the corresponding Icr errors is ctr = 0.019. 



(23) 



IV. RESULTS 



A. CPL parameterization 



We will present the results of CPL parameterization in this subsection. It is found that^^^.^ = 461.254, 
465.440, and 466.100 for SNIa, SNIa+BAO, and SNIa+BAO+CMB data, respectively. Moreover, we 
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FIG. 1: Reconstructed w(z) and f(z) from Snia data alone (upper panel) and SnIa-i-BAO-i-CMB data(lower panel) 
using the CPL ansatz. Red lines stand for the best-fit values and green lines stand for the Icr CL. 

reconstructe w{z) and /(z) by using CPL parameterization, as is shown in figure [T] The upper panel of 
figure \T\ shows reconstructed w{z) and /(z) from SnIa data alone. From this panel, it is seen that w and 
/ (i.e. DE density p\) will rapidly decrease along with the increase of redshift z, and the behavior of DE 
significantly deviate from the cosmological constant A at Icr CL. It should be mentioned that our results are 

n 

consistent with some other works. For example. Figure 2 of lilTH shows that the best-fit value of wi is rather 
negative. This means that the EOS w of DE will also be quite negative at high-redshift region. 

The lower panel of figure [T] shows reconstructed w(z) and /(z) from SnIa-i-BAO-i-CMB data. For this 
case, the result is quite different, and the ACDM model is still consistent with observational data at Icr CL. 
The difference from these two panels of figure [T] shows the existence of the tension among different types 
of observational data, which was recently discussed in [32]. Their result of the CPL parametrization shows 
that for the Union compilation when CMB and BAO data is added there is a Ax^ = 1.063, while for the 



Constitution set A^^ = 4.846, which is much larger. Their result is consistent with ours. Besides, in 13311 
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the authors get (wq, wi) - (-0.95 0.41 ;^Jj™) for the CPL parametrization from SnIa+BAO data, which 
shows the existence of tension between Snia and BAO data compared with the upper pannel of figure [T] 

As mention above, for fitting SNIa data alone, we find that both w and pa will rapidly decrease along 
with the increase of redshift z by using CPL parameterization. This result implies that the behavior of 
dark energy might be very different in different slices of redshifts, and inspires us to separate redshifts into 
several bins and to investigate the piecewise constant w model and the piecewise constant p\ model. 



B. ACDM2 and XCDM2 Model 

Let us discuss the cases in which redshifts are separated into 2bins. Hereafter we will call 2 bins piece- 
wise constant w model as XCDM2 model, and will call 2 bins piecewise constant pA model as ACDM2 
model. 

Figure|2]shows the. x^min versus redshift z. for XCDM2 and ACDM2 models, where only the Constitution 
SNIa sample is used in the analysis. It is seen that the curves of x^„in of these two models are very similar, 
and both these two models achieve their minimal at the discontinuity points of redshift zv = 0.331. We 
will analyze the XCDM2 and ACDM2 model more explicitly at this discontinuity point. The result is as 
follows. 

For XCDM2 model, using the SNIa data and the discontinuity point z\ = 0.331, we find that the best-fit 
value and Corresponding Icr CL of the model parameters are Q„,o = 0.466 ;J;° g22' ~ ~l-118^Q 4jg, and wi 
with the best-fit value -508.170 and the upper Icr CL value -15.2. For this model, = 456.574 (As a 
comparison, for the CPL parametrization to the same data, theXmin ~ 461.254. The XCDM2 model lower 
the Xmin 4-32). Notice that the EOS of the second redshift slice wi <c -1, this means the DE density pA 
should decay very quickly at the range z > 0.331. 

For ACDM2 model, using the SNIa data and the discontinuity point zi = 0.331, we find that the best- 
fit value and Corresponding Icr CL of the model parameters are Q^o = 0.461^Q gg2 and q = (5.407 x 
10^)^0462' corresponding X^n^ - 456.893. Notice that the best-fit Q,„o and the, Xmin "^^^y close to the 
result of the XCDM2 model. We also plot the Icr and 2cr error bars for the ACDM2 model in figure [3l It 
is seen that there is a deviation from cosmological constant A over Icr confidence level, which is similar to 
the results of CPL parameterization and XCDM2 model. 

To further verify this conclusion, we construct three toy models to mimic the step function, as is listed in 
table in and analyze them using the SNIa data. It is found that, for toy model 1, the best-fit model parameters 
are n„,o = 0.457 and ^ = 3.734x 10"^' , corresponding to^^ . = 457.805; for toy model 2, the best-fit model 
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FIG. 2: versus redshift z for XCDM2 model and for ACDM2 model. Only the Constitution SNIa sample is used 
in the analysis. 
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FIG. 3: Estimates of Pa/pao of the ACDM2 model. Only the Constitution SNIa sample is used in the analysis. This 
figure shows a deviation from the ACDM over \<j CL. 



TABLE I: Three toy models that mimic the step function 



Model 


Toyl 


Toy2 


Toy3 


/(z) 






1 






l+gk-Jo)/'- 




457.805 


457.499 


456.771 
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FIG. 4: The evolution of DE density pA of three toy models. 

parameters are Q„,o - 0.462, ^ = 2.644 x 10"^' and s - 3.108, corresponding to xf„i„ - 457.499; for toy 
model 3, the best-fit model parameters are Q,„o = 0.462, z.o = 0.329, and v = 4.514 x 10~'^ corresponding 
to^^.^ = 456.771. Notice that the numbers of free parameters of these three toy models are less than or as 
same as that of CPL parameterization, so both these three models perform better than CPL parameterization 
in fitting the Constitution set. we also plot the evolution of f{z) of these three models in figure |4j where the 
best fit parameters of each model are adopted. It is seen that all these three models have a sharp decrease in 
Pa at z ~ 0.33, which is consistent with the ACDM2 model and the XCDM2 model. 

As is shown in figure [51 we also take the influence of BAO and CMB data into account. Using the 
SNIa+BAO data, the best-fit ACDM2 model has a discontinuity point zi = 0.975, and its best-fit value and 
Corresponding Icr CL of the model parameters are n„,o - O.m'^'^ '^^^ and q = -^^^-t^^T ' corresponding 
to xl,j„ = 461.906. Using the SNIa-i-BAO-FCMB data, the best-fit parameters of ACDM2 model has a 
discontinuity point zi = 0.859, and its best-fit value and Corresponding Icr CL of the model parameters 
are Q^o = 0.283;|;qq23 ^'^'^ ^2 - ^■^'^'^lojs2' corresponding to xlii„ - 463.552. Therefore, although the 
step function model in which pA rapidly decreases at redshift z ~ 0.331 is very powerful in fitting the SNIa 
data, it is not the best model to fit the combination of SNIa, BAO and CMB data. This shows the existence 
of the tension among different types of observational data |32]. This result is similar to the result of CPL 
parameterization. 

For a complete comparison, using SNIa, SNIa-i-BAO, and SNIa-i-BAO-i-CMB data, we Ust x^jn, 
xiiJdof, and BIC for XCDM2 model and ACDM2 model in table [III table [Till and table |lVl respec- 
tively. It is seen that the xl,i„ of ACDM2 model are very close to those of XCDM2 model. However, since 
the number of free parameters of ACDM2 model is one fewer than that of XCDM2 model, the xLml^^f 
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FIG. 5: xlnn versus redshift z for ACDM2 model, where SNIa, SNIa+BAO, and SNIa+BAO+CMB data are used, 
respectively. 

TABLE II: The x\,i„ for XCDM2 model, ACDM2 model and CPL parameterization 



Model 


SNIa 


SNIa+BAO 


SNIa+BAO+CMB 


XCDM2 


456.574 


462.500 


464.775 


ACDM2 


456.893 


46L906 


463.552 


CPL 


46L254 


465.440 


466.100 



and BIC of ACDM2 model are smaller than those of XCDM2 model. This means the ACDM2 model 
performs better than the XCDM2 model in fitting the observational data. To make a further comparison, 
we also list the )^„iiJdof, and BIC for CPL parameterization, respectively. One can see that, for 
SNIa, SNIa+BAO, and SNIa+BAO+CMB data, the . and BIC of ACDM2 model are 4.361, 3.534, and 
2.549 smaller than those of CPL parameterization, x\,iJd.of of ACDM2 model are 0.011, 0.009, and 
0.006 smaller than those of CPL parameterization. Therefore, the ACDM2 model presents a significant 
improvement for the same data set obtained using the CPL ansatz. 



C. ACDM3 and XCDM3 Model 

Then, let us turn to the cases in which redshifts are separated into 3 bins. In table |Vl we list ;^f^ and 
BIC for 3 bins piecewise constant w (XCDM3) model and 3 bins piecewise constant pA (ACDM3) model. 
Comparing with the cases of 2 redshift bins, the reduction of . is very small, but the increase of BIC is 
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TABLE III: The Xm Jdof forXCDM2 model, ACDM2 model and CPL parameterization 



Model 


SNIa 


SNIa+BAO 


SNIa+BAO+CMB 


XCDM2 


1.162 


1.174 


1.177 


ACDM2 


1.160 


1.169 


1.171 


CPL 


1.171 


1.178 


1.177 



TABLE IV: The BIC for XCDM2 model, ACDM2 model and CPL parameterization 



Model 


SNIa 


SNIa+BAO 


SNIa+BAO+CMB 


XCDM2 


480.510 


486.446 


488.731 


ACDM2 


474.845 


479.865 


481.519 


CPL 


479.206 


483.399 


484.068 



quite obvious. Again, one can see that the Xmin ACDM3 model are very close to those of XCDM3 model, 
but the BIC of ACDM3 model are much less than that of XCDM3 model. Therefore, for same number of 
bins, a piecewise constant pa model always performs better than a piecewise constant w model. This shows 
the advantage of using pA, instead of w, to probe the variation of DE. 

We will not discuss the result of the ACDM3 and XCDM3 model explicitly. Here we only mention that 
the second bins of the ACDM3 and XCDM3 models also show a sharp decrease in p\, as is consistent with 
the CPL parameterization, the LCDM2 model and the XCDM2 model. We get the minimal = 455.112 
at two discontinuity points z. = 0.331 and z = 0.975, and plot best-fit values and Icr CL of pa/pao for the 
ACDM2 model in figure [6l As is consistent with figure [T] , figure |4] and figure [3l this figure also shows a 
visible decrease of pA at z > 0.331 and a deviation from the ACDM in Icr CL. Notice that the Error Bar of 
the third bin is too big to get any convincing information, the reason is that there are only 20 supernovaes 
in this bin. 



TABLE V: Thexli„ (BIC) for XCDM3 model and ACDM3 model 



Model 


SNIa 


SNIa+BAO 


SNIa+BAO+CMB 


XCDM3 


455.713 (491.617) 


459.598 (495.517) 


462.028 (497.962) 


ACDM3 


455.112 (485.032) 


460.109 (490.041) 


461.059 (491.004) 
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FIG. 6: Estimates of pa/paq of the ACDM3 model. Only the Constitution SNIa sample is used in the analysis. 

V. SUMMARY 



In this work, we explore the cosmological consequences of the recently released Constitution sample 
of 397 SNIa. By revisiting the CPL parameterization, we find that, for fitting the Constitution set alone, 
the behavior of DE significantly deviate from the cosmological constant A, where w and pA will rapidly 
decrease along with the increase of redshift z. Inspired by this clue, we separate the redshifts into different 
bins, and discuss the models of a constant w or a constant pA in each bin, respectively. It is found that for 
fitting the Constitution set alone, w and pA will also rapidly decrease along with the increase of z, which is 
consistent with the result of CPL model. Moreover, a step function model in which pA rapidly decreases 
at redshift z ~ 0.331 presents a significant improvement (Ax^ = -4.361) over the CPL parameterization, 
and performs better than other DE models. We also construct three toy models that mimic this process, and 
find a DE model that close to the step function model will be favored by the Constitution SNIa data. By 
plotting the error bars of the piecewise constant DE density of ACDM2 and ACDM3 model, we show that 
there is a deviation from cosmological constant A at about Icr confidence level; this may arise from some 
biasing systematic errors in the handling of SNIa data, or more interestingly from the nature of DE itself. 
In addition, for models with same number of redshift bins, a piecewise constant pA model always performs 
better than a piecewise constant w model; this shows the advantage of using pA, instead of w, to probe the 
variation of DE. 
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